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Abstract—New procedures are developed in this article for estimating pa- 
rameters of the Pearsonian Type I curve which are particularly adaptive 
to factors influencing the pattern of age-specific fertility rates. It is shown 
that with this model the number of parameters required for the graduation 
and simulation of these rates can be reduced to only three—total fertility 
rate, mean and modal ages of fertility. The reduction in the number of 
fertility parameters offers considerable operational and analytical ad- 
vantages, and makes the Pearsonian Type I curve particularly appropriate 
for the construction of a parametric model for fertility projections. In light 
of the results of empirical tests based on fertility data for Canada, the 
model’s potential for birth projections appears quite promising. 


INTRODUCTION 


Statisticians and demographers have 
_ shown considerable interest in the Pear- 
sonian Type I curve, primarily as a 
means of graduation of age-specific fer- 
tility rates (Wicksell, 1931; Mitra, 1967; 
Keyfitz, 1968; Avery, 1970). Although a 
few authors have also recognized its 
potential for birth projections (Tekse, 
1967; Stone, 1970), to the best of our 
knowledge no one has yet systematically 
investigated this aspect. This article has 
two objectives: first, to simplify the 
procedures for estimating the parameters 
of the equation of the Pearsonian Type I 
curve, which is employed in the gradua- 
tion of age-specific fertility rates; and 
second, to investigate the feasibility of 
utilizing this curve as a foundation in 
the construction of a parametric model 


for birth projections. The two objectives 


are closely related; the first objective is 
in fact a prerequisite of the second. 
However, a model which is satisfactory 
for graduation purposes might be un- 


satisfactory for projections. For example, 
higher order moments may help to per- 
fect the graduation, but they are often 
inappropriate for forecasting because 
they are difficult to interpret in demo- 
graphic terms. One of the high-priority 
goals of demographers is to develop a 
projection model of increased operational 
and analytical capability. However, if 
this goal is to be attained, it is essential 
that a procedure be developed which 
would rely on a limited number of demo- 
graphically meaningful fertility param- 
eters, and simultaneously, one which 
would permit age-specific fertility rates 
to be satisfactorily reproduced. The 
derivation of a procedure which would 
satisfy these two criteria is the ultimate 
goal of this article. 

In the first section, the Pearsonian 
Type I curve is briefly presented, and a 
justification is provided for its selection 
as a mean of reproducing age-specific 
fertility rates. In the second section, 
there is a short review of the procedures 
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usually used for deriving constants in 
the equation of the Pearsonian Type I 
curve, followed by a presentation of new 
methods which have been developed ta 
estimate the constants. These new pro- 
cedures are tested in the third section by 
using Canadian historical data on fer- 
tility, and in the final section, attempts 
are made to explore the ways in which 
the Type I function could be utilized for 
birth projections. 


PRESENTATION AND JUSTIFICATION 
OF THE TYPE I CuRVE 


With the origin at the mode, the equa- 
tion of the Type I curve is: 


Y= yl t/a)” () —.2/an)" 4 


where y, is the modal ordinate and 
—a<x< a. 
Also, 


m,/d, = M/a. (2) 


In formula (1), y, is the modal fertility 
rate, which in human populations falls 
between 20 and 30 years of age; a, and 
@ jointly determine the reproductive 
interval; and m, and m, determine the 
shape of the fertility curve. When m, 
and my, are approximately equal, the 
curve approaches a normal distribution, 
but it is positively skewed when m, < 
ms. Modern reproductive patterns of 
population are characterized by the lat- 
ter skewed distribution; this is demon- 
strated in Appendix Table A-1, which 
presents the parameters of a fertility 
distribution obtained from data for 
Canada. One can see, as a first indication 
of skewness, that the modal age is con- 
sistently smaller than the mean age, and 
that the measure of skewness, 8), is posi- 
tive throughout the 1926-1969 pcriod. 
However, the selection of a particular 
Pearsonian curve to describe the fertility 
experience of any given country cannot 
be based merely on impressionistic views 
and judgments based on the parameters 
mentioned above. Indeed, K. Pearson 
(Elderton, 1930) has argued that the 
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criterion for this selection should be 
based on the calculation of k: 


k= Bi (Be ae 3)’/[4 (2B. ai 6) 


(4B, a 38:)]. 


Pearson recommends that the Type I 
curve should be used when k has a nega- 
tive value. 

It can be seen in Appendix Table A-1 
that for Canada k has a negative value, 
and consequently, that the selection of 
this: curve is justified here. It should be 
noted, however, that k assumes progres- 
sively higher negative values. In view 
of these shifts in the age pattern of 
Canadian fertility, it would be worth- 
while trying some of the other curves in 
the Pearsonian system of curves, such 
as the “normal curve” or its proxy, 
Type II, to represent the age pattern 
of fertility in earlier years and the tran- 
sitional Type III curve to represent the 
age pattern of fertility in more recent 
years. Alternatively, outside of models 
based on the so-called gamma distribu- 
tion, only the Mazur-type model (Mazur, 
1963) has proved to be appropriate in 
the case of relatively high and less 
asymmetrical fertility (Tekse, 1967). 

A comparative study of fitting various 
curves is outside the scope of this paper. 
Instead, we shall concentrate our atten- 
tion on the curve which, based on the 
criterion of k, appears to be most suit- 
able for the period under observation and 
which, by implication, will probably con- 
tinue to conform to Canada’s fertility 
patterns in the near future. Tests con- 
ducted by Mitra (1967) and Avery 
(1970), using fertility materials for 
nations throughout the world, have re- 
vealed that the Pearsonian Type I func- 
tion is among the most suitable mathe- 
matical functions developed so far for 
approximating observed age schedules of 
fertility. Unfortunately, in view of data 
constraints, their fitting of the curve had 
to be restricted to fertility rates by five- 
year age groups, and consequently a de- 
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tailed analysis of the goodness of fit was 
not possible. 

Our next task is to discuss the various 
methods by which the constants in for- 
mula (1) can be derived. 


MeEruHops oF DERIVING CONSTANTS 


The method most commonly used for 
fitting formula (1) to the actual age- 
specific fertility rates is the method of 
moments. However, the number of mo- 
ments used depends on the particular 
procedure followed. In the Elderton 
(1930) procedure, the constants a, de, 
m,, and mz in equation (1) are calcu- 
lated from the first four moments of the 
frequency distribution, that is, from the 
mean, variance, skewness, and kurtosis. 
Later, Mitra (1967) developed a pro- 
cedure that reduces the required number 
of moments to the first two, but this 
procedure assumes a fixed age interval 
of fertility. Thus, with the fertility in- 
terval set at ages 15 through 50, 


a; se ag = 35. (3) 


With equations (2) and (3), the required 
number of independent parameters is re- 
duced to only two, and the solutions are 
obtained from the following equations: 


m, + mM, = pa’ (ai 1S My’) / Ma a: 


(4) 
and 
m, = [(m, + m2+ 2)/(ai + a2)Jm.’ — I, 
(5) 


where p,;’ represents the mean counted 
with the starting point at age 15 and 
wy, the variance. Equations (2) and (3) 
can then be used to determine a, and ay. 
As for y., the following can be applied: 


yo = N/[(a, + a2)B(m, + 1, m, + 1)] 


*[m,™"»m,™")/[(m, + 9 Reales (6) 


In equation (6), N represents the sum 
of the age-specific fertility rates and B 
represents the Beta function. In fact, y, 
may be regarded as a multiplier that 
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equalizes the sum totals of the observed 
and graduated distributions. 

The above procedure is based on 
equating the first two moments of the 
observed distribution with those of the 
theoretical distribution. Alternatively, 
the procedure for calculating the con- 
stants may be further simplified by using 
the first moment of the distribution, the 
mean, and some other measure, such as 
the mode. A conventional method may be 
utilized to estimate the mode, or one 
could simply substitute the midpoint of 
the single-year age interval with the 
highest age-specific fertility rate. In this 
case, given the value of the mode or a; 
with the origin at the start of the curve, 
the remaining parameters can be ob- 


tained from the following: 


(a; 'a,) —.a), (7) 
Mz = A2(a, + Az — 2y’)/[(a, + ae) 
“Gi 1d), (8) 


a2 


and 


mM, = (a;/d2)Mp. (9) 


Although the reduction in the number 
of moments may entail some loss in 
goodness of fit, it presents some definite 
analytical and operational advantages. 
These advantages will be examined in 
greater detail later, but at the present 
time let us merely emphasize that the 
model’s parameters should respond in a 
manner that is consistent with the trend 
in fertility rates. In the absence of this 
consistency, it would not be possible to 
relate the variations in the distributions 
to the variations in the parameters. Con- 
sequently, the model would be of little 
practical value. Given the multitude of 
possible parameters, it is extremely dif- 
ficult to relate, for example, the effect of 
a reduction in ‘the total fertility rate 
on the modal age, the fertile age span, 
etc. The restrictions on the fertility in- 
terval in the second procedure and on 
the fertility interval and modal age in 
the third procedure make the interpre- 
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tation of the remaining parameters rela- 
tively simple and perhaps more imean- 
ingful. 

Other procedures for deriving or per- 
fecting the estimation of the constants 
may be suggested. One such procedure, 
based on the fertile age range, the modal 
age, and the modal fertility rate, also 
seems reasonable and worthy of in- 
vestigation. With the origin at the start 
of the curve, and using the relationship 


m,/d, = M2/A2 


= (m, + m2)/(a, + a) = C, (10) 
equation (6) for the modal fertility rate 
can be rewritten as: 


N 


(a) (ey an 


where N for the single-year distribution 
is equivalent to the total fertility rate. 
The constant C, obtained through an 
iterative procedure, can then be used in 
(10) to estimate the remaining param- 
eters. 

A different approach that could be 
used to improve the fit has been sug- 
gested by Keyfitz (1968), namely, itera- 
tion techniques based on a_ descent 
method. This approach proved useful in 
fitting the Gompertz curve to Canadian 
data on cumulative fertility (Murphy 
and Nagnur, 1972), and with due adap- 
tation it could be applied to fertility 
data in the form of age-specific distribu- 
tions. 

The two latter procedures are men- 
tioned only to indicate further develop- 
ments to improve the fit. At this point 
we were not able to fully investigate 
their potential for projections. Instead 
we have concentrated on the first three 
procedures and have tested them against 
empirical data, the results of which are 
to be presented in the next section. 

For simplicity in future references, 
these three procedures will henceforth be 
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referred to by the following codes shown 
in parentheses: 


a) Method based on the first 
four moments 

b) Method based on the first 
two moments 

c) Method based on the mean 
and the mode 


(4M) 
(2M) 


(1M). 


The numerals 4, 2, and 1 stand for the 
number of moments on which a par- 
ticular method is based. 


EMPIRICAL VERIFICATIONS OF THE 
PROCEDURES AND THE MopDEL 


Whatever may be the theoretical justi- 
fication of a particular method, any 
decision concerning its suitability for 
practical application must rest on some 
empirical test, such as the comparison 
of derived estimates with actual figures. 
Before proceeding with the testing, a 
brief discussion of the data base is in 
order. 

These tests are performed using data 
on fertility by single years of age for 
Canada available since 1926. Except for 
the personal affiliation of one of the 
authors, there is no particular reason 
for selecting this country. Canadian data 
on fertility are assumed to be fairly 
good, but not completely free from re- 
porting errors. This is probably true, par- 
ticularly for the earlier years which, 
among other things, are characterized 
by gross irregularities in the modal age 
group. There seems to be no theoretical 
argument in favor of the oscillations ob- 
served at alternate ages in this age 
group. The extent to which the remain- 
ing ages have been affected by this kind 
of error is not known, and no attempt 
has been made to adjust the data for 
such biases. However, the magnitude and 
direction of these biases are not expected 
to greatly influence the characteristics 
of the distributions on which the esti- 
mates of the parameters are based. 

Only period fertility data are used 
here. Data by cohorts are reconstructions 
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from the period data, supplemented by 
extrapolations from truncated cohorts. 
While these cohort data are useful for 
general types of analysis, as such they 
do not warrant enough confidence to 
base on them the assessment of the 
validity of the procedures being tested 
here. 

For convenience of space, the years 
1926, 1931, 1941, 1951, 1961, and 1969 
have been selected for statistical com- 
parison. During this time period there 
were considerable variations in both the 
level and age pattern of fertility. In this 
paper, estimates for each year within the 
1926-1969 period are provided only for 
method (4M) (Appendix, Table A-2). A 
priori, one may expect the closeness of 
fit to improve as the number of moments 
used is increased. It follows that method 
(4M) sets some standard against which 
estimates obtained by other methods 
can be gauged. 
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We shall compare first the estimated 
fertility measures and parameters with 
the actual ones, and then examine the 
goodness of fit for the overall fertility 
age schedule, obtained by the three pro- 
cedures described in the preceding sec- 
tions. 


Estimated and Actual Fertility Measures 


Of all the parameters of the Type I 
distribution, the one that is most mean- 
ingful in this context seems to be the 
modal age. Since method (2M) is re- 
stricted by a fixed fertility interval, the 
estimate of the modal age will not be 
unique and will depend on the specific 
choice of that interval. The results based 
on two alternative fertility intervals of 
15-50 and 17-50 years by (2M) along 
with those calculated by (4M) are com- 
pared with observed modal ages in Table 
1. It should be noted that the observed 
values presented in the table were ad- 


TaBLeE 1—Hstimates of Modal Age for Canada Using Methods (4M), (2M), and (1M): 1926, 
1931, 1941, 1951, 1961, and 1969 


Modal Age Obtained from Method 


Total 
CSSA ON ca (1M) Fertility 
Year (4M) (15-50) (17-50) (observed) (per 1000) 
(1) (2) (3) (4) (5) (6) 
1926 28.1 b 28.9 Zhe 28.0° 3356 
(17.0-48.5) 
293k 27.0 b 28.5 27.6 ZT? 3201 
(16.9-49.3) 
1941 2053 b 27.6 26.4 2752 2824 
(17.2-49.6) 
1951 25.5 b 26.6 25.4 Fhe P| 3480 
(17.0-49.9) 
1961 24.2 b 23e4 24.3 24.0 3857 
(17.2-50.4) 
1969 24.0 : “ts pd» oe il ys 2410 


(16.6-53.0) 


a-Fixed fertility interval. 


b-Estimate of the fertile age range based on (4M). 
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justed for irregularities which appear to 
have been caused by random fluctuations 
and age misreporting. These adjustments 
were made by tracing a freehand curve 
through the points of observed modal 
ages. It should also be recalled that 
method (1M) has no built-in mechanism 
for the calculation of modal age, and 
that the latter has to be given in order 
to enable us to use this method for 
deriving other parameters. Methods (4M) 
and (2M), with an age interval of 17-50, 
estimate the modal age fairly well, 
whereas method (2M), with a fertility 
range of 15-50, tends to overestimate it 
by about one year. 

The figures in parentheses in column 
(2) are the estimates of the fertile age 
range based on (4M), and it can be seen 
that the lowest (16.6) and the highest 
(53.0) limits occur in 1969. The lower 
limit is no larger than 17.2 years, and 
the smallest value of the upper limit is 
48.5 years. The latter appears to vary 
more than the former. In general, the 
derived measure of reproductive period 
is not quite in accord with the actual one 
observed among human _ populations. 
Such inconsistencies are unavoidable 
when all the parameters are estimated 
from the moments of distribution. 

As for the modal fertility rate shown 
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in Table 2, no clear pattern emerges 
from the differences between the esti- 
mates obtained by various methods. 
They all appear to estimate modal fer- 
tility fairly well, with the exception of 
(1M) (15-20), which tends to under- 
state it considerably. 

Parameters m, and mz, which were 
derived by various procedures, are shown 
in Tables 3 and 4. As could be expected, 
the values of these parameters are, and 
should be, independent of the total fre- 
quency since their derivation is depen- 
dent primarily upon the distribution of 
the relative frequencies. Accordingly, 
changes in the parametric values are re- 
flections of the changes in the pattern of 
the distributions themselves. Parameter 
m, appears to be particularly sensitive 
to such changes. This is reflected in the 
high coefficients of correlation shown in 
Table 5, which indicates the association 
between m, and various parameters of 
the age pattern of fertility, such as mean 
and modal ages, variance, skewness, and 
kurtosis. Over the period of observation, 
m- has exhibited a rather steady upward 
trend, which merely seems to reflect the 
downward shift in the age pattern of 
fertility in Canada. During the same 
period, when one disregards relatively 


Taste 2.—Estimates of Modal Fertility Rate for Canada Using Methods (4M), (2M), and 
(1M): 1926, 1931, 1941, 1951, 1961 and 1969 


Modal Fertility Rate Obtained from Method Ob 
served 
Modal 
oa $$) gereiiity 
Year (4M) (15-50)* 17-50)" (15-50) (17-50)* Rate 
1926 178 181 181 152 184 188 
1931 174 175 175 150 181 183 
1941 159 159 160 147 173 168 
19351 205 203 205 172 204 210 
1961 239 234 236 195 231 257 
1969 157 150 155 130 p Leh 168 


a-Fixed Fertility interval. 
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Taste 3.—Estimates of m: for Canada Using Methods (4M), (2M), and (1M): 1926, 1921, 
1941, 1951, 1961, and 1969 


m4 Obtained from Method 


(2M) (1M) 
Year (4M) (15-50)* (17-50)? (15-50)* (17-50)? 
1926 90 io) 98 <o3 1505 
1931 395 LSD 97 89 1.09 
1941 3 1.49 89 lee ke Lazo 
1951 “od 1.45 82 ota 81 
1961 he bray 73 .64 64 
1969 95 Lin39 72 «16 74 
a-Fixed fertility interval. 
minor fluctuations (see Appendix, Table Goodness of Fit 


A-2), m; has remained relatively stable; 
apparently, it has been less influenced by 
the changing age pattern of fertility. 
However, its apparent high (negative) 
correlation with the total fertility rate, 
shown in Table 5, is difficult to explain. 

In light of the results presented in 
Tables 3 and 4, it is also quite apparent 
that m, and m, must be very sensitive to 
the estimation procedure. The estimates 
obtained by the methods differ among 
themselves and in relative terms; the 
magnitude of these differences are con- 
siderable. 


So far, attention has been focused on 
the agreement between the observed and 
estimated fertility parameters. We shall 
now see to what extent the derived fer- 
tility age distribution agrees with that of 
the observed. 

The classical approach of measuring 
goodness of fit is to calculate X?, the 
magnitude of which depends in part upon 
the relative difference between the ob- 
served and graduated values. This 
method was considered inappropriate for 
a situation characterized by a large num- 
ber of intervals, for which some of the 


Taste 4—Estimates of m2. for Canada Using Methods (4M), (2M), and (1M): 1926, 1931, 
1941, 1951, 1961 and 1969 


m9 Obtained from Method 


(2M) (1M) 
Year (4M) (15-50)° (17-50)* (15-50) (17-50)* 
1926 1.67 25 2.00 1.40 aba 
1931 1.93 2.46 2.08 hed 2.28 
1941 2.12 2.68 225 2.14 2.82 
1951 2.43 2.92 2.43 1.81 2.40 
1961 2.64 342 2257 1.86 2.39 
1969 3.71 3.42 2.80 aids 2.84 


a-Fixed fertility interval. 
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(4M) estimates were zero, namely, those 
at the beginning and the end of the fer- 
tility curve. By comparing the observed 
and expected distributions, one can ob- 
tain overall index values in a number of 
ways. One such index, the index of dis- 
similarity (A), is obtained by reducing 
the two distributions to percentage form 
and summing only the positive differ- 
ences between corresponding percentages. 
A given value of this index indicates the 
percentage of observations which must be 
redistributed among intervals in order 
that the two distributions become iden- 
tical. It is easy to see that A can range 
from zero to 100 and that its magnitude 
depends on the choice of class intervals. 

In Table 6, the overall error appears 
to be relatively small. As could be ex- 
pected, the index values are consistently 
smallest for the (4M) method, but con- 
sidering the simplicity and logical con- 
sistency of the other methods, the values 
are not very high. Hence, in general the 
model seems to be quite satisfactory, at 
least as a first approximation of a gradu- 
ation formula for fertility rates. Only 
method (1M), with a fertility age in- 
terval of 15 to 50, results in more con- 
sequential error. 

In order to obtain some indication of 
the kind of age-associated biases gen- 
erated by the model, the deviations of 
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the derived from the observed age-spe- 
cific fertility rates have been calculated 
and depicted in Figure 1. The outstand- 
ing feature of this figure is the fact that 
for all observed years there is a very 
definite pattern in the age biases. The 
model tends to alternate between under- 
estimating and overestimating the fer- 
tility rates for a series of successive age 
categories. Similar cyclical patterns of 
age biases are generated by the other 
two procedures. 

If the model is to be used for projec- 
tions, it is important to ascertain the 
structural conditions of fertility for 
which the model performs best. To this 
effect, correlation coefficients between the 
index of dissimilarity and the measures 
of the age pattern of fertility were cal- 
culated. As Table 5 reveals, all of the 
correlation coefficients are significant, 
but the highest coefficient (.70) is the 
one obtained with criterion k, which can 
be viewed as an index of the combined 
effect of skewness and kurtosis on the 
fertility curve’s overall shape. It follows 
that as k increases, that is, as its nega- 
tive values become increasingly larger, 
the error generated by the model be- 
comes progressively larger. Consequently, 
after k reaches a certain level, an alter- 
native curve such as the transitional 
Type III curve mentioned earlier may 


Taste 6—A Values for Canada Using Methods (4M), (2M), and (1M): 1926, 1931, 1941, 
1951, 1961 and 1969 


vis Values Obtained from Method 


2M) (1M) 
Year (4M) (15-50)4 (17-50)? (15-50)o (Pi=50)2 
1926 2.76 4.08 a7 7 8.94 2.80 
1931 2.66 3.78 2.59 715 o37 
1941 2.32 4.30 2.39 5133 4.51 
1951 2.28 4.03 pe 9.42 2.41 
1961 2.97 5.38 3.30 11.34 4.03 
1969 3.48 5,23 3.99 10.39 4.13 


a-Fixed fertility interval 
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provide a better fit than the Type I 
curve. However, the negative values for k 
which were obtained from Canadian fer- 
tility data are not nearly large enough 
to justify the utilization of the Type III 
curve. The errors in the results which are 
generated by each of the three methods 
differ from each other in magnitude, but 
not in pattern. Therefore, although the 
coefficients discussed here refer to the 
results derived by method (4M), the in- 
ferences to which they give rise apply to 
each of the methods. 


POTENTIALS FOR BIRTH PROJECTIONS 


In the preceding sections, the Type I 
curve has been examined as a device for 
the graduation or reproduction of ob- 
served age-specific fertility rates. In the 
present section, an attempt will be made 
to extend its application to the field of 
projections. We shall first demonstrate 
how this can be done, and then illustrate 
how the operational and analytical ad- 
vantages of the model based on this 
curve make it a more effective device 
for projecting birth series than the con- 
ventional methods. 

In trying to adapt the Type I curve to 
projection requirements, it is helpful to 
distinguish between parameters which 
refer to the level of fertility and those 
which refer to the age distribution of fer- 
tility. It should be recalled that in 
formula (1) the former are represented 
by the modal fertility rate, y,, and the 
latter by the parameters aj, d2, m,, and 
Mo. 

The fertility level may be projected in 
_one of two ways. One could project the 
modal fertility rate, in which the model 
in its formulation by Pearson is stated, 
or alternatively, one could project the 
total fertility rate. Although there is a 
very high correlation (.93) between the 
modal and total fertility rates, on sta- 
tistical and analytical grounds it seems 
preferable to project the latter rates. 
This is because the historical series of 
modal fertility are influenced by irregu- 
larities which are probably attributable 
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to age misreporting and random fluctua- 
tions. Moreover, the utilization of total 
fertility rates permits one to conduct 
analysis in much more depth than is pos- 
sible with the modal fertility rate. For 
example, the cohort approach can be 
employed to explain variations in the 
total fertility rate by referring to changes 
in family size and shifts in the birth 
timing of successive cohorts of women 
(Ryder, 1969). Furthermore, future fer- 
tility trends can be forecast by analyzing 
parity distribution data and information 
obtained from surveys on intended family 
size (Siegel and Akers, 1969). These ad- 
vantages of employing the total fertility 
rate as the basis for projecting future 
fertility levels cannot be attained if one 
attempts to make projections in terms 
of modal fertility; in this latter case the 
forecaster can rely only on past trends 
as a guide to the future. 

Let us now turn our attention to the 
projection of fertility age distributions. 
In this case one can directly project 
either the dependent parameters, a, a2, 
Mm,, and ms, or the independent param- 
eters, that is, the moments and other 
measures which are associated with fer- 
tility frequency distributions. The selec- 
tion of one of the three retained proce- 
dures as the basis for calculating the 
dependent parameters determines the 
choice of the appropriate independent 
parameter. 

In our view, there are primarily two 
prerequisites which must be met if a 
particular set of parameters is to be 
selected for incorporation into a model 
for projections, and these are analytical 
suitability and amenability to a demo- 
graphically meaningful interpretation. 
Thus, a model which yields a “good fit” 
but does not simultaneously meet the 
prerequisites mentioned here may be 
satisfactory for graduation purposes, but 
it would be totally inadequate for pro- 
jection purposes. According to these 
criteria, dependent parameters which 
provide mathematical descriptions but 
do not provide meaningful demographic 
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interpretations do not qualify for se- 
lection, At this point the difficulties we 
have encountered earlier in this arti- 
cle in offering a demographically mean- 
ingful interpretation of the behavior of 
m, and mz should be recalled. To a cer- 
tain extent, moments of higher orders 
suffer a similar fate. 

Among the three listed procedures, the 
one which, in our opinion, best satisfies 
the criteria of “rationality” and demo- 
graphic meaningfulness of parameters is 
method (1M), which requires the knowl- 
edge of merely two simple and relatively 
easily understood fertility measures, 
namely, mean and modal ages of fer- 
tility. Moreover, since these two mea- 
sures are highly correlated (.98), for all 
practical purposes one may confine the 
analysis and projection to only one of 
them. Whereas it is true that procedure 
(1M) does not provide quite as close a 
fit as do other procedures such as (4M), 
the resulting loss in precision is more 
than compensated for by the increased 
operational and analytical advantages 
which are possible with the former 
method. 

The accuracy of a projection based on 
any particular method can only be de- 
termined ex post facto. However, for 
the sake of argument, let us momentarily 
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assume that all the independent param- 
eters had been correctly projected for 
the 1926-1970 period, and that each of 
the models developed in this paper had 
been used to reproduce the total annual 
numbers of births during this period. 
Given these assumptions, any differences 
between the actual totals of annual 
births and those derived from the models 
must necessarily be due to the models. 
The ratios of the number of births esti- 
mated by the different models to the ac- 
tual numbers of births for Canada in 
selected years are shown in Table 7. One 
can see that the ratios are all close to 
unity. This is true even for method 
(1M), which relies on only three fertility 
measures, namely, the total fertility rate 
and the mean and modal ages of fertility. 
Consequently, in the calculation of the 
total annual number of births in the 
indicated period, each of the three 
models performs almost to perfection. 


CONCLUSION 


Additional research is needed before a 
parametric model of fertility projections 
of the type outlined here can be made 
fully operational. This paper contains, 
nevertheless, important components for 
such a model. It is shown that fertility 
distribution by age can be derived by 


TaBLe 7.—Ratio of Estimated to Actual Number of Annual Births for Canada Using Methods 
(4M), (2M), and (1M): 1926, 1931, 1941, 1951, 1961, and 1969 


Ratio of Estimated to Actual Number 
of Births Obtained from Method 


RRO ORT Cs. re pull LA a GS Ptr aiada 
Year (4M) (15-50)* (17-50)* (15-50)3 (17-50)? 
1926 1.001 986 .999 1.005 1.000 
1931 1.001 1.001 1.000 1.010 .998 
1941 .998 .997 .998 1.000 .996 
1951 .998 .997 997 984 .996 
1961 1.003 1.005 1.001 1.013 1.002 
1969 1.008 1.006 1.011 1.018 1.009 


a-Fixed fertility interval 
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mathematical functions from only a _ that cannot be achieved with conven- 
limited number of parameters that need tional procedures making direct use of 
to be projected. The paper offers alter- age-specific fertility rates in projecting 
native ways of calculating these param- births. Furthermore, the model is a 
eters; it tests their results against fer- powerful labor-saving device, for in de- 
tility data for Canada; and it points to fining fertility in terms of a mathemati- 
potentials for further developments inthe cal function, computers can be utilized 
areas of fertility projections. Reduction to perform many of the involved opera- 
of the age schedule of fertility to only a tions. 

few meaningful parameters makes pos- 

sible an in-depth analysis to a degree APPENDIX 


Taste A-1—Parameters of Canadian Fertility Distribution by Age, 1926 to 1969 


Total Measure Measure 
Fertility Observed Mean Std. of of 
Rate Modal Modal Age of Devia- Skewness kurtosis Criterion 

Year (Per 1000) Fertility Age@ Fertility tion By Bo K 
1926 3360 188 27.98 SORnS 6.61 .059 2,285 -0.029 
1927 3320 181 27.94 s0nLi 6.62 (059 2.284 -0.029 
1928 3300 196 27.86 30.04 6.62 062 2.286 -0.030 
1929 3220 176 PME SH 29.90 6.60 077 eo le -0.037 
1930 3280 182 PM Ages) 29.87 6.58 082 2,521 -0.040 
1931 3200 183 27.70 29.86 6.54 083 2.348 -0.042 
1932 3090 b7fs3 27.40) 29.96 6.57 .975 Te SVAN -0.038 
1933 2870 163 IML ASN) 29.98 6.56 063 2.340 -0.032 
1934 2800 154 Pai Eres yee 30.08 6.54 054 Yara yee -0.028 
1935 2750 158 PREY 30.00 6.54 .057 2.345 -0.030 
1936 2700 156 2742 29.96 6.52 -058 20337 -0.029 
1937 2640 148 PEL ARNG 29.83 6251 069 2.344 -0.035 
1938 2700 by Pat 3X0} 29.67 Ooo 084 2351 -0.043 
1939 2650 ihSyl 27 fe eas) 29.61 6.50 09) 259 -0.046 
1940 2760 159 27.19 29.36 6.47 <5 2.388 -0.058 
194) 2820 168 Zan 29.15 6.42 143 Behe. -0.072 
1942 2950 17) 26.94 292 6.39 eee: 2.446 -0,073 
1943 3000 183 26.80 29.15 6.34 M35 2.454 -0.071 
1944 3000 AS 26.60 29.30 6.37 114 2.415 -0.059 
1945 3000 — 177 26.42 2938 6.40 vive 2.394 -0.056 
1946 3360 209 Aloe) 29.01 6.33 i LSS 2.464 -0.080 
1947 3580 225 26.05 28.72 6.29 -184 2.500 -0.095 
1948 3420 at! 2ST 28.67 Ge .181 2.498 -0.094 
1949 3440 213 25470 28.60 6.29 181 aeols -0.095 
1950 3430 213 25.50 28.64 6.36 -109 Zee -0.100 
1951 3480 210 2533 28.45 6.28 3287 Drone -0.099 
1952 3620 222 25615 28.35 6.26 197 2.536 -0.104 
1953 3700 225 24.96 28.29 6.25 e £202 e535 -0.106 - 
1954 3810 242 24.78 28.25 6.27 avaalte 2.548 -0.112 
1955 3820 245 24.60 28in22 6a25 220 2551 -0.114 
1956 3850 254 24.40 Poh Wd 6.23 w2o2 2.584 -0.123 
1957 3930 253 246.23 28.00 6625 e242 2.587 -,126 
1958 3880 247 24,14 27399 One 249 2005 -O.131 
1959 3950 Zi 24.10 27.87 6.19 . 263 2.629 -0,140 
1960 3910 255 24.06 Zoo 6.17 270 2.636 -0.143 
1961 3860 257 24602 earl ThE 6.16 . 284 22655 -0.151 
1962 3770 205 24.00 20eS (arin 291 2.669 -0.155 
1963 3690 249 Zoe OT PAs jee fis) 6.10 . 286 ZO -0.155 
1964 3520 236 23.94 27.80 6.1) 284 2.68) -0.156 
1965 3160 Zu2 23.90 27.76 6.14 276 2.668 -0.151 
1966 2820 193 23.86 PAT Moy 4 (Ae ils} 296 2.700 -0.163 
1967 2590 182 23.84 Daal 6.07 mee Phe Lsy4 -9.182 
1968 2440 170 23.80 27' S28 5.98 341 2O37: -0.209 
1969 2410 168 PL Saha) Blie2e 5.98 341 2.837 -0.208 


-a~Freehand smoothed values of observed modal] age of fertility rate. 
Source: Data for the calculation of this table are taken from the annual reports of Vital 
Statistics, 1926-1969, Canada, Dominion Bureau of Statistics. 
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Taste A-2.—Parameters Estimated by Elderton’s Method: 1926-1969 


a 
ee DeROSIER SEE SES AES PERRO EEE ee ee oi ea 


Modal y 

Year Age. 0 ay ee tl yi) 

FORK isk wah. SESS oA Ue ee ee 
1926 28.08 177.90 11.10 20.50 0.90 1.67 
1927 28.04 175.33 pie v2 20.55 0.90 1.66 
1928 Soyee » 174.35 10.99 20.64 0.89 1.66 
1929 27.59 171.83 10.59 21.26 0.87 1.75 
1930 27.51 175.95 10.43 21.41 0.86 et 
1931 27.60 173.79 10.68 21.68 0.95 1.93 
1932 27.83 166.81 11.03 valet 0.99 1.93 
1933 28.09 161.01 11.61 21.20 1.09 1.99 
1934 28.34 151.25 11.93 20.82 is 1.97 
1935 28.26 149.00 11.99 21.08 aw, 2.05 
1936 28.21 146.61 11.88 20,89 1.14 2.01 
1937 27.86 143.78 11.30 21.24 1.06 1.99 
1938 27.43 147.35 10.64 21.66 0.96 1.95 
1939 27.24 144.93 10.34 21.79 0.92 1.93 
1940 26.73 153.20 9.69 22.48 0.86 2.00 
1941 26.30 159.39 9.14 23.27 0.83 yen: 
1942 26.38 167.70 9.34 7337 0.89 2522 
1943 26.58 year 9.62 23.25 0.96 2.93 
1944 26.87 169.89 9.97 D253 0.97 2.19 
1945 26.84 168.40 9.88 b2- 37 0.92 2.07 
1946 26.17 194.03 8.98 23.50 0.86 2.25 
1947 25.62 209.85 8.34 24.09 0.79 2.29 
1948 25.59 199.17 8.46 24.17 0.80 2.30 
1949 25.63 201.82 8.65 24.30 0.86 244 
1550 25.55 199.72 8.55 24.74 0.83 2,01 
1951 25.45 205.07 8.54 24.43 0.85 2.43 
1952 25.28 214.82 8.34 24.59 0.83 2.45 
1953 25.14 219.89 8.15 24.55 0.80 2.41 
1954 24.99 226.43 7.96 24.90 nF7 2.41 
1955 24.88 228.13 7.76 24.88 0.74 2.38 
1956 24.81 232.15 7.76 25.34 0.77 a, 53 
1957 24.55 236.81 7st 25.45 0.73 2.47 
1958 24.51 236.10 7.49 25.61 0.75 2255 
1959 24.39 241.98 7.35 25.96 0.74 2.62 
1960 24.29 241.05 7.19 25.95 car 7,64, 
1961 24.17 239.29 7.01 26.24 0.70 2.64 
1962 24.15 PAG 27 6.91 26.26 0.71 2.68 
1963 2495 231.65 7.07 26,27 0.74 2.75 
1964 24.38 220.58 707 26.49 0.78 2.85 
1965 24.36 196.56 7.38 26.40 0.79 2.8% 
1966 24.14 176.60 7.20 26.94 0.77 2.90 
1967 23.91 165.30 TaGh 27.54 0.79 3.09 
1968 24.04 159.44 7.39 28.98 0.95 3.71 


Pearsonian Type | Curve 
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